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PREDICTION  OF  HELICOPTER  ROTOR  POWNWASH 
IN  HOVER  AND  VERTICAL  FLIGHT 


SUMMARY 

~A  method  is  presented  for  the  calculation  of  rotor  downwash  in  hovering  and  vertical 
flight.  The  blade  is  represented  by  a lifting  line.  The  rotor  wake  is  simulated  by  a vortex 
sheet  and  a series  of  rolled-up  root  and  tip  vortices,  similar  in  form  to  that  of  a classical 
fixed-wing.  The  concept  of  rectangularization  of  the  rotor  wake  is  used  in  obtaining  a 
formula  for  the  normal  component  of  induced  velocity.  Wake  contraction  based  on  experi- 
mental data  is  introduced  into  the  calculations.  Numerical  calculations  have  been  performed 
for  two  rotor  configurations,  viz.  S- 58, and  a model  rotor.  Computed  induced  velocities  and 
blade  loadings  are  compared  with  the  available  flight  data.  The  comparison  shows  the 
present  simple  method  yields  satisfactory  results. 
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NOTATION 

a lift-curve  slope 

A wake  contraction  constant 

b number  of  blades 

r chord  of  the  blade  at  station  r 

C ratio  0-5  ac,'(/?2  - Ri) 

Cl  section  lift  coefficient 

Ct  rotor  thrust  coefficient 

h distance  between  successive  vortices 

h,  ratio  2sh  (fl_>  Ri) 

dl  lift  on  an  element  of  the  blade  dr 

m radial  station  along  the  blade 

M total  number  of  radial  stations  along  the  blade 

N number  of  rolled-up  vortices 

r radial  coordinate  along  the  blade 

Rr  position  of  root  vortex 

Rt  positions  of  tip  vortex 

Ri  inner  radius  of  the  blade 

R;  outer  radius  of  the  blade 

.s  number  identifying  the  particular  vortex  under  consideration 

dT  thrust  on  an  element  of  the  blade  dr 

l',  vertical  velocity  ( • ve  up) 

»•(  induced  velocity  ( i ve  down) 

v0  mean  induced  velocity 

v induced  velocity  at  P due  to  vortex  sheet 

v,  induced  velocity  at  P due  to  sth  rolled-up  vortex 

H'  incident  velocity  (see  eqn.  5) 

H'r  weight  of  the  aircraft 

a incidence  angle 

r circulation  of  station  r 

To  circulation  round  the  centre  of  the  aerofoil 

Tj  integration  variable  (see  eqn.  18) 

f*  spanwise  coordinate  (see  eqn.  16) 

d pitch  of  the  blade  at  station  r 
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1.  INTRODUCTION 

Rotor  performance.  blade  motion  stability.  blade  stresses,  vibrators  hub  forces,  aircraft 
stability  and  control,  and  aircraft  flight  simulation  are  some  of  the  areas  which  require  a know- 
ledge of  the  non-unttorm  downwash  at  the  rotor  In  the  (light  simulation  studies  which  are  our 
main  concern  the  rotor  downwash  is  required  in  the  calculation  of  rotor  thrust  and  other  rotor 
forces  and  moments  It  is  also  needed  in  the  calculation  of  aerody namic  forces  on  the  helicopter 
fuselage,  attached  weapons  and  any  other  externally  slung  birdies.  Probably  one  of  the  most  difficult 
problems  facing  the  rotary  wing  aerody  namicist  at  the  present  lime  is  the  determination  of  this 
induced  velocity  or  rotor  downwash  As  each  new  generation  of  helicopters  becomes  more 
sophisticated,  increased  emphasis  is  put  on  improving  the  rotor  induced  velocity  calculations. 

I he  assumption  that  the  distribution  of  induced  velocity  over  the  disc  is  uniform  is  thought 
to  be  satisfactory  in  forward  Might  when  the  Mow  through  the  rotor  is  large.  In  forward  flight, 
with  forward  tilt  ot  the  rotor  disc,  the  blades  are  moving  away  from  the  wake,  therefore  feeling 
less  of  its  influence*.  In  the  hovering  case,  however,  the  wake  beneath  the  rotor  may  be  expected 
to  play  a prominent  part  in  determining  the  induced  vehvcity.  In  this  report,  this  important 
Might  condition  is  the  one  analysed  A simple  but  rigorous  analytical  technique  which  will  aid 
in  understanding  and  interpreting  the  complex  relationships  between  rotor  circulation,  down- 
wash  and  blade  motion  will  be  ol  great  value  in  helicopter  flight  simulation  studies.  This  problem 
has  been  the  subject  of  considerable  work  over  many  years,  l'he  calculation  techniques  have 
evolved  from  simple  momentum  theory  based  on  actuator  disc  theory,  through  the  classical 
blade  element  momentum  method,-  to  wake  modelling  pnvccdures.  The  advent  of  the  high 
speed  electronic  digital  computer  has  made  possible  the  straightforward  approach  of  tracing 
the  vortex  filaments  trailed  by  each  blade,  and  integrating  the  Biot-Savart  relation  to  obtain 
induced  velocities  The  wake  modelling  methods  progressed  to  the  sophisticated  and  compu- 
tationally demanding  free  wake  analysis,*1  7 then  retreated  to  the  simpler  prescribed  wake11  13 
methods.  In  the  prescribed  wake  method,  the  wake  geometry  is  specified,  on  the  basis  of  experi- 
mental data,  as  a function  of  rotor  configuration  and  thrust  level,  through  simple  analytical 
expressions.  In  the  free  wake  analysis  the  mutual  interaction  of  the  circulation  distribution, 
downwash  and  blade  motion  is  included  Many  theoretical  models  based  on  the  above  free 
wake  concept  have  been  developed  in  I urope.  United  States  of  America  and  United  Kingdom. 
T ue  majority  of  these  models  rely  exclusively  on  discrete  vortex  filament  methods  which  are 
based  on  extensive  digital  computation.  These  techniques  are  necessary  when  detailed  rotor 
blade  airloads,  vibration  analysis,  and  acoustic  properties14  are  required  However,  because 
of  the  extensive  numerical  computations  involved,  they  oiler  only  limited  physical  insight,  and 
lack  the  versatility  of  simpler  methods  l ven  when  using  the  most  sophisticated  models,  con- 
vergence in  the  hovering  range  is  very  doubtful3. 

The  aim  of  the  following  analysis  is  to  develop  a vortex-wake  method  for  determining 
the  induced  vehvcity  of  a hovering  helicopter  The  method  uses  a simple,  realistic  wake  model 
based  on  experimental  data  and  Willmcr's,-S"1  mathematical  simplification  to  avoid  lengthy, 
complex  computations  The  basic  geometry  of  the  vortex-wake  model  is  shown  in  Tigure  I. 
In  this  model  each  blade  is  represented  bv  a lilting  line  The  vorticity  trailed  by  a rotor  blade 
leaves  the  trailing  edge  in  a sheet,  which  rolls  up  as  a pair  of  concentrated  tip  and  root  vortices. 
In  order  to  reduce  the  computational  tune.  Willmer's  principle  of  rectangulari/ation  is  applied 
to  the  above  wake  model  Using  this  vortex-wake  model,  the  calculated  downwash  distribution 
is  compared  with  available  theoretical  and  experimental  results 

In  developing  the  theory,  consideration  has  been  limited  to  aircraft  having  a single  rotor. 
No  restriction  has  been  placed  on  the  number  of  blades  which  the  rotor  may  have,  however. 
The  aircraft  has  been  assumed  to  be  in  steady  vertical  Might  or  hover.  The  fuselage  interference 
on  the  rotor  wake  is  not  incorporated  It  is  also  assumed  that  the  Muid  is  inviscid  and  incom- 
pressible. The  simplicalions  employed  in  formulating  the  mathematical  nu*del  restrict  its 
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application  to  points  in  the  rotor  plane,  but  the  model  could  he  extended  to  include  points  in 
other  regions  of  interest  The  calculated  rotor  downwash  ma>  be  used  as  a first  approximation 
in  (he  estimation  of  aerodynamic  forces  on  the  helicopter  fuselage  as  well  as  any  weapons 
attached  to  it. 


2.  ROTOR  WAKE  VORTKTTY  DISTRIBUTION 


The  underlying  aerodynamic  theory  required  will  be  better  understood  if  the  equivalent 
fixed  wing  problem  shown  in  f igure  2 is  first  examined.  Consider  a wing  advancing  at  constant 
angle  of  attack  in  a uniform  airstream  By  virtue  of  its  geometrical  angle  of  attack,  the  aerofoil 
will  generate  circulation  and  consequently  lift  Assuming  for  simplicity  that  this  circulation  is 
constant  along  the  span,  it  must  leave  the  wing  at  the  tips  and  trail  downstream  since,  by  one 
of  the  fundamental  laws  of  hydrodynamics,  circulation  cannot  end  abruptly,  but  must  continue 
back  in  a closed  circuit  to  the  starting  vortex  generated  at  the  beginning  of  the  motion.  After 
the  motion  has  continued  for  a sufiicieni  length  of  time,  this  starting  vortex  may  be  assumed 
to  be  at  an  infinite  distance  I t orn  the  aerofoil  and  its  effects  may  be  neglected.  In  practice  of 
course  the  circulation  is  not  constant  along  the  span  and  therefore  the  trailing  vortices  are 
distributed  in  a sheet  all  along  the  span  Consider  the  case  w here  the  geometrical  angle  of  attack 
changes  as  the  wing  advances  It  is  evident  that  the  circulation  will  change  but.  since  the  total 
circulation  must  remain  constant  in  the  case  of  the  ideal  fluid,  there  must  be  a counter  vortex 
in  the  wake  corresponding  to  the  change  in  circulation  of  the  aerofoil  This  wake  vort icily  is 
generally  referred  to  as  the  shed  vorticity  to  differentiate  it  from  the  trailing  vorticity  which 
occurs  both  in  the  steady  state  condition  and  when  the  angle  of  attack  of  blade  is  changing 
In  steady  level  flight,  the  variations  in  the  angle  of  attack  are  small  so  that  the  magnitude  of 
the  resulting  shed  vorticity  will  also  be  small.  Its  contribution  to  the  induced  velocity  is  neglected. 
It  is  also  well  established  that  the  vortex  sheet  behind  an  aerofoil  is  unstable  and  cannot  persist.17 
The  sheet  tends  to  roll  up  somewhat  like  a sheet  of  paper  as  shown  in  figure  2 Thus,  at  a sufficient 
distance  behind  the  aerofoil,  a section  of  the  wake  in  a plane  perpendicular  to  the  direction  of 
motion  would  show  two  cylindrical  vortices  whose  distance  apart  is  less  than  the  span  (refer 
to  Fig  4) 

In  Figure  I the  corresponding  problem  for  the  rotor  is  presented  Obviously  the  trailing 
wake  does  not  disappear  to  infinity  but  returns  underneath  the  rotor  in  a spiral  form"'  and 
the  very  nature  of  this  spiral  prevents  the  development  of  simple  solutions  for  the  downwash 
which  it  induces  at  the  blade  Because  of  this  complexity,  generalized  solutions  will  probably 
never  be  possible  for  the  problem  However,  numerical  solutions  for  specific  cases  are  possible 
utilizing  high  speed  digital  computers  These  solutions  are  frequently  sufficient  for  engineering 
applications 

Figure  I indicates  pictorially  the  complexities  of  attempting  to  obtain  a complete  represen- 
tation of  the  rotor  wake  To  solve  this  difficult  problem  numerically.  Willmer  made  two 
simplifying  assumptions  The  first  is  that  conditions  change  sufficiently  slowly  to  allow  a quasi- 
static approach.  The  second  is  that  the  spiral  wake  is  amenable  to  rectangulari/ation.  Consider 
the  trailing  wake  from  a blade  (Fig  ')  Willmer  argued  that  the  radius  of  curvature  of  the  vortex 
sheet  and  rolled  up  vortices  is  large  enough  (especially  Tor  the  important  outer  parts  of  the 
wake)  for  the  wake  elements  (vortex  sheet  and  rolled  up  vortices)  to  be  regarded  as  straight. 
He  also  assumed  that  only  those  parts  of  the  wake  which  arc  near  to  the  reference  blade  are 
important  Consistent  with  this  idea  one  may  allow  the  rows  of  vorticity  to  extend  to  infinity 
in  order  to  achieve  mathematical  simplification  Thus  the  wake  attached  to  the  reference  blade 
is  assumed  to  be  a straight  sheet  extending  back  to  infinity,  while  those  of  other  blades  are 
assumed  to  be  doubly  infinite  rolled-up  vortices  as  shown  in  Figure  4.  The  positions  of  these 
vortices  depend  on  the  mean  flow  velocity  through  the  rotor  disc,  and  the  number  and  relative 
positions  of  the  blades  that  shed  them  Once  the  numbers  and  positions  of  the  wake  elements 
had  been  chosen,  the  induced  velocity  at  the  reference  point  could  be  calculated  by  an  extension 
of  Glauert's  wing  theory  Positioning  of  these  vortices  forms  one  of  the  important  aspects  in 
the  calculations  of  induced  velocity.  The  axial  distance  h between  the  vortices  is  calculated  using 
the  average  induced  velocity  through  the  rotor.  A method  based  on  the  experimental  data  is 
used  to  determine  the  radial  locations  This  is  discussed  in  the  following  section. 
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3.  WAKE  GEOMETRY 


The  most  important  factor  in  the  accurate  calculation  of  rotor  downwash  today  is  the  use 
of  realistic  wake  geometry  Failure  of  most  of  the  earlier  classical  performance  prediction 
methods  was  attributed  to  the  use  of  an  uncontracted  wake.  In  the  following  pages  we  briefly 
describe  the  wake  geometry  used  in  the  present  analysis. 

Concurrently  with  the  rotor  inflow  methods  based  on  a theoretical  wake  geometry,  other 
methods  have  been  developed  based  on  empirical  wake  models  Most  of  these  methods  are 
directed  toward  the  rotor  hover  condition  in  that  this  condition  is  the  one  most  influenced  by 
wake  distortion  effects.  The  requirement  for  a method  employing  an  accurate  prescribed  wake 
model  derived  from  experimental  wake  data  was  concluded  by  Jenney,  et  a/11  who  found  that 
the  rapid  contraction  of  the  slipstream  under  a hovering  rotor  places  the  vortex  system  sufficiently 
close  to  the  rotor  blades  that  it  causes  significant  changes  in  the  distribution  of  induced  velocity. 
In  the  mid  I950's,  Gray*  developed  a semi-empirical  method  for  the  wake  of  a single  bladed 
rotor  based  on  experimental  wake  geometry  data  obtained  from  smoke  visualization  tests.  More 
recently.  Landgrebc.®  and  Kocurek  and  Tangier10  have  conducted  a series  of  smoke  tests  which 
confirm  Gray's  results  and  give  more  details  of  the  wake  geometry.  Kocurek  and  Tangier 
measured  wake  shapes  for  26  model  rotor  configurations  described  in  Table  I. 


TABLE  I 

Model  Rotor  Configurations1  u 


Untwisted  Rotors 


Blade  linear 

Aspect 

mm 

Radius 

Chord 

ii  R, 

twist  (deg) 

ratio 

R»  (m) 

c (mm) 

(m  s) 

0 

71 

0 407 

57-2 

152  5 

0 

13  5 

1.  2.  3 

0 554 

41  4 

152-5 

0 

18  2 

1.  2.  3.4 

0 680 

37-3 

152-5 

Twisted  Rotors 

no 

13  5 

1.  2.4 

— 

41  4 

1 52  5 

113 

8 0 

1.  2.4 

1 

57-2 

152  5 

16  0 

71 

1,  2.4 

1 

57-2 

152-5 

18  0 

80 

1.  2,  4 

0-458 

57-2 

152-5 

24  0 

71 

1.  2.4 

0-407 

57-2 

152-5 

All  the  blades  had  a rectangular  plan  form  and  a NACA  001 2 airfoil  section. 
Landgrebc  obtained  the  wake  geometry  characteristics  for  all  combinations 
of  the  rotor  parameters  given  in  Table  2. 
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tabu:  2 

Values  of  Primary  Tes!  Parameters* 


Number  of  blades 

2,  4,  6.  8 

Blade  linear  twist  (deg) 

0,  8,  16 

Blade  as|>ect  ratio 

15  6,  18  2 

Rotor  tip  speed  (in  s) 

160,  185,  214 

Collective  pitch* 

0 to  maxt 

* Value  at  r equal  to  0 75  R: 
t l>etermined  by  operating  \tall  limits. 


The  above  two  How  visualization  studies  revealed  that  the  wake  radial  contraction  occurs 
in  an  exponential  manner  with  increasing  wake  azimuth.  ^ It  is  characterized  by  a rate  parameter, 
k,  and  an  effective  minimum  nondimensional  radius.  4.  The  generalized  wake  radial  coordinate 
is  given  by 

' 4 ■ (I  4 ) e (i) 

Aj 

For  the  stable  near-wake  region,  the  radial  coordinates  may  be  determined  by  substitution  of 
the  following  combination  of  4 and  A in  the  above  equation  ( i >. 

4 0 78 

ft)  145  27  ()(';•  l.andgrebe 

A < (ii) 

[d  ll  (('/•)*  Kocurek  and  Tangier 

Landgrebe’s  linear  variation  of  A agrees  with  the  non-linear  representation  of  Kocurek  and  Tangier 
in  the  medium  Cr  range.  However,  in  the  case  of  the  latter  authors  the  contraction  rate  tends 
towards  zero  as  the  blade  is  unloaded.  From  the  above  equation  (u)  it  is  concluded  that  for  a 
given  azimuth,  the  radial  coordinate  of  the  wake  is  a function  of  disc  loading  (or  Cr)  and 
decreases  with  increasing  loading,  l or  a lived  disc  loading  the  radial  coordinates  are  independent 
of  tip  speed,  number  of  blades,  blade  twist,  and  aspect  ratio.  Taking  a mean  Cr  value  of  0 005 
and  h'R equal  to  0-058,  the  wake  geometry  is  plotted  in  Figure  5.  The  Figure  shows  how  rapidly 
the  wake  contracts  under  the  rotor:  the  contraction  is  practically  complete  within  a distance  of 
20  to  40  percent  of  the  rotor  radius. 

Using  the  above  described  wake  flow  model  and  the  wake  contraction  rate  given  by  Kocurek 
and  Tangier,  a mathematical  analysis  is  developed  in  the  following  section.  The  variation  of 
the  circulation  with  radius  and  the  resulting  down  wash  distribution  is  represented  by  a Fourier 
scries.  Unfortunately  , since  the  loading  of  a rotor  blade  is  very  dilferent  from  the  elliptic  loading 
of  a conventional  fixed  wing,  a large  number  of  spanwise  points  are  needed  for  acceptable 
accuracy.  This  means  a large  number  of  terms  in  the  Fourier  series,  and  this  makes  for  lengthy- 
computing  times. 

4.  MATHEMATICAL  FORMULATION 

When  a wing  experiences  lift,  (he  pressure  on  its  upper  surface  must  be  less  than  on  its 
lower  surface.  It  follows  from  Bernoulli’s  theorem  that  the  velocity  is  higher  over  the  upper 
surface  than  over  the  lower  surface,  and  we  should  therefore  expect  that  the  circulation  round 
a circuit  surrounding  the  aerofoil  would  be  non-zero.  Following  the  lifting-line  practice,  the  lift 
of  an  element  of  the  blade  dr  and  at  a distance  r from  the  origin  is17 

dl-  Fp W dr  (I) 

When  the  inflow  angle  (the  direction  of  flow  relative  to  the  blade)  as  shown  in  Figure  6 is  small 
we  can  write 

dT  rPH  dr  (2) 
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In  the  above  relationship  neither  thrust  nor  induced  velocity  arc  known.  To  obtain  a further 
relationship,  we  now  consider  the  lilt  characteristics  of  the  blade  where  the  blade  is  regarded 
as  an  aerofoil.  The  calculations  follow  closely  the  standard  methods  of  aerofoil  theory  but  the 
rotor  analysis  is  simplified  considerably  because  the  blade  incidence  and  inflow  angles  are 
usually  sufficiently  small  that  the  familiar  small  angle  approximations  may  be  made. 

Consider  an  element  of  blade  of  chord  r with  width  dr  at  a radius  r from  the  axis  of  rotation. 
The  geometric  pitch  angle  of  the  blade  element  relative  to  the  plane  of  rotation  is  6,  the  climbing 
speed  is  l/r  and  the  local  induced  velocity  is  »•,.  The  direction  of  the  flow  relative  to  the  blade 
makes  an  angle  $ (usually  called  inflow  angle)  with  the  plane  of  rotation,  and  is  given  by  (refer 
to  Fig.  6) 


tan 

(IV  • Vi);Or 

(3a) 

The  lift  on 

or  <P  s 

the  blade  element  is 

(L’c  ■ vi)!Or  for  small  <t> 

(3b) 

dL 

} P W '1  c Cl  dr 

(4) 

where 

W 

(((/,.  • v,)-  ( 

(5a) 

and 

H * 

Qr 

(5b) 

for  small  <P 

Let  us  suppose  that  the  lift  slope,  a.  of  the  section  is  constant  so  that,  if  the  section  incidence  a 
is  measured  from  the  no-lift  line,  we  can  write 

Ci.  a*  - - <P)  (6) 

The  elementary  lift  is  now 

dL  \p  fi'-r4  a(8  — <t>)  c dr  (7) 

Using  eqn.  (3b)  and  assuming  Jl.  5 JT 

dT  J p uc  Q2r-  (*?  ((  V • v()  Qr]  dr  (8) 

Thus  from  eqns.  (2)  and  (8)  the  relation  between  the  circulation  and  the  induced  velocity  becomes 

F J ac  ( BQr  t V n)  (9) 


As  in  the  fixed-wing  case,  the  circulation  round  the  rotor  blade  will  vary  across  the  span.  Between 
the  points  and  (>j  • <h))  of  the  span  of  the  blade,  the  circulation  changes  by  the  amount 
— drjdr)  <h)  and  hence  a trailing  vortex  of  this  strength  springs  from  the  element  d^  of  the  span 
as  shown  in  Figure  7. 

The  normal  induced  velocity  at  any  point  P along  the  span  must  be  obtained  as  the  sum  of 
the  effects  of  all  the  trailing  vortices  immediately  behind  the  blade.  The  induced  velocity  due 
to  the  trailing  vortex  sheet  is’® 


V sherl 


I 

4tt 


r2  dr 


<b) 


r - v 


(10) 


ti 


The  next  contribution  comes  from  the  rolled-up  vortices  stacked  below  the  rotor  blade 
near  the  root  and  tip  as  shown  in  Figure  8.  The  induced  velocity  perpendicular  to  the  blade 
at  P due  to  an  clement  of  rolled-up  vorticity  in  the  s'th  wake  at  Q and  Q'  is’® 


I (r  R,)ra 

2, r(r  R,)i  • s*h* 

so  that  the  total  for  the  whole  wake  is 
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Thus  ihe  total  induced  velocity  at  P (serpendicular  to  the  tip-path  plane,  which  is  given  by 
that  due  to  the  immediate  wake  plus  the  sum  of  the  contributions  of  the  rolled-up  vortices 
beneath  the  rotor,  is 


lt  ' htrt  — t* 

From  eqns  (9).  (10),  (12)  and  (13)  the  relation  for  the  circulation  /'becomes 


(13) 
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(14) 

We  now  express  the  circulation  in  Fourier  senes  form: 

/’  iJR.{R.  R>)  ^ Am  sin  mt>  (15) 

m 1 

where  l>  is  defined  by  the  relation 

I R:  ■ R,\  iR-  R\\ 

' ( , ) ( 2 ) co>  ' ’ 

The  integral  in  eqn.  ( 14)  is  singular  in  nature  To  evaluate  it  we  use  the  well  known  result17 


cos  >n<i>  M n sin  mO 
cos  <b  cos  0 sin  0 


(17) 


J 0 

The  relation  between  the  variable  y in  eqn.  (14)  and  <t>  in  eqn.  (17)  is 


The  next  step  is  to  determine  I\i,  Rt  and  Rr  We  use  the  theory  developed  by  Glauert  for 
a fixed  wing  aerofoil  shown  in  Figure  9.  The  trailing  vortex  sheet  behind  an  aerofoil  is  unstable 
and  rolls  up  into  a pair  of  vortices  whose  distance  apart  /'  is  rather  less  than  the  span  / of  the 
aerofoil.  There  is  no  established  method  the  calculate  to  strength  of  the  trailing  vortices  springing 
from  an  unsymmetrically  loaded  blade.  Here  we  assume  that  the  strength  of  these  trailing  vortices 
will  be  equal  to  the  magnitude  of  the  circulation  r « round  the  centre  of  the  aerofoil.  If  the  form 
of  load  distribution  across  the  span  is  represented  as  in  equation  (15)  by  the  Fourier  Sine  series, 
then 
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Using  eqns.  (15)  to  (18),  eqn.  (14)  reduces  to 
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From  the  wake  geometry  we  also  know  the  value  Rt.  Using  the  relation  (21)  we  could 
calculate  Rr  the  position  of  the  rolled-up  root  vortex.  Substituting  eqn.  (20),  eqn.  (19)  reduces  to 
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Assuming  the  vertical  velocity  of  rolled-up  vortices  is  equal  to  the  average  induced  velocity 
plus  the  climb  velocity,  the  vertical  distance  between  successive  vortices  is  given  by 
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(29) 


Selecting  finite  values  for  m and  .v  say  M and  A respectively,  equation  (22)  yields  a set  of 
linear  simultaneous  algebric  equations  from  which  one  can  evaluate  the  Fourier  coefficients. 
After  this,  calculation  of  induced  velocity  and  thrust  coefficient  is  a simple  matter.  From 
equations  (9),  (15),  (16),  (24)  and  (28)  the  induced  velocity  is  given  by 
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and  from  equations  (I),  (15)  and  (5b)  the  blade  loading  distribution  is 
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(31) 


Integrating  the  above  equation  along  the  blade,  multiplying  by  the  number  of  blades  and  dividing 
by  pA&Rj  we  get  the  thrust  coefficient  of  the  rotor 


CT  - 


(i  - m r i )-  f 
4 [2 


1 -f 


At 


(32) 


On  substituting  A equal  to  zero  in  the  above  analysis  we  obtain  the  results  for  the  hovering  case. 
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5.  NUMERICAL  EXAMPLES 


An  iterative  process  is  used  to  calculate  induced  velocity  and  thrust.  An  arbitrcv  downwash 
distribution  and  thrust  coefficient  are  first  assumed.  The  wake  geometry  is  generate*  tsing  these 
initial  values.  Then  we  formulate  a set  of  simultaneous  algebraic  equations  with  the  Fourier 
coefficients  as  unknowns.  These  equations  are  solved  by  using  the  standard  numerical  algorithm 
based  on  C rout's  factorisation  method.  From  known  Fourier  coefficients,  new  values  of  induced 
velocity  and  thrust  coefficient  are  calculated  thus  completing  the  first  iterative  cycle.  The  process, 
which  is  shown  in  Figure  10,  is  repeated  until  convergence  of  v«  and  Ct- 

Calculations  were  performed  to  determine  the  induced  velocity  of  two  helicopter  rotors. 
The  aircraft  is  assumed  to  be  in  steady  hovering  flight.  The  important  parameters  for  these 
rotors  are  given  in  Table  3. 


TABLE  3 

Important  Parameters  for  the  Two  Rotors 


Rotor 

' 

Radius  (m) 

1 

Chord 

1 

Tip 

! 

No.  of 

! 

Blade 

Linear 

Aero- 



(mm) 

Speed 

Blades, 

Root 

Twist 

Foil 

Ri 

R> 

c 

(r/min) 

b 

Pitch 

Rate 

Section 

S-582" 

1 -37 

8-50 

222 

4 

13-9° 

-8  0° 

NACA 

0012 

Model  Rotor11 

0-088 

0-58 

j 51-0  1 

3 

150° 

—8-0° 

NACA 

0012 

The  weight  of  the  model  aircraft  was  about  6-8  kg.  The  weight  of  the  S-58  aircraft  varied 
between  5085  and  5360  kg.  In  the  numerical  calculations  M and  N were  set  to  20  and  12  respec-  f-. 

lively.11  In  both  cases,  lift  slope  is  taken  to  be  5.73.  From  published  experimental  results910  Ct  £ 

was  initially  assumed  to  be  0 005.  An  initial  value  of  v©  was  calculated  using  the  momentum  , 

theory  equation2: 

*«.  •}[  U-  1 (U?  i 2 H-V/mt /?!)•]  (33)  , 

With  the  above  input  data,  the  iterative  process  is  commenced.  Calculated  values  of  v©  are 
plotted  as  a function  of  iteration  number  in  Figure  1 1.  As  we  can  see  from  this  figure,  the  solution 
converges  very  rapidly. 


6.  COMPARISON  WITH  FLIGHT  DATA  AND  WITH  OTHER  THEORIES 


The  results  shown  in  Figure  12  indicate  that  the  incorporation  of  wake  contraction  signifi- 
cantly alters  the  radial  distribution  of  induced  velocity.  This  is  particularly  true  in  the  blade 
tip  region  where  the  rolled-up  vortex  passes  close  to  the  blade.  Comparison  of  calculated  values 
with  the  experimental  data  is  limited  by  the  amount  of  data  available.  In  fact  there  are  no 
experimental  data  on  induced  velocity  in  the  rotor  plane  at  the  present  time.  Hence  computed  blade 
loading  is  compared  with  measured  blade  loading20  in  Figure  13.  Here  once  again  we  seethe 
effect  of  wake  contraction  on  the  results.  One  should  be  cautious  in  comparing  the  present 
mathematical  model  with  other  theories,  as  all  theories  necessitate  many  simplifying  assumptions. 
Therefore  any  comparison  would  rest  upon  the  significance  of  the  simplifying  assumptions  and 
hi.w  they  affect  the  calculated  values.  Figure  15  shows  a comparison  with  Jenney,  Olson  and 
Landgrebe’s  sophisticated  mathematical  model.  The  results  compare  well  over  the  greater  part 
of  the  blade  span  except  near  80",,  of  the  blade  radius.  This  could  possibly  be  due  to  the  simpli- 
fication of  the  wake  representation  by  planar  rather  than  curved  surfaces. 

Results  presented  in  Figures  12-14  also  compare  the  present  theory  with  Willmer's  theory 
and  flight  data.  Included  in  these  figures  are  the  calculation  both  with  and  without  wake  con- 
traction. It  can  be  seen  from  Figure  13  that  the  present  method  more  closely  predicts  the  blade 
loading  in  the  tip  region  where  the  loading  is  very  high.  This  is  mainly  due  to  the  incorporation 
of  wake  contraction. 
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7.  CONCLUDING  REMARKS 


The  method  described  provides  a rapid,  computer-oriented,  numerical  technique  for 
evaluation  of  rotor  downwash.  Results  presented  for  two  different  rotors  indicate  that  the 
variation  of  induced  velocity  along  the  radius  of  the  rotor  is  highly  non-linear.  Correlation  of 
the  computed  and  measured  results  is  quite  good,  considering  the  simple  wake  geometry  used. 


1 
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FIG.  4.  TRAILING  VORTEX  SHEET  AND  ROLLED  UP  VORTICES  NEAR  BLADE  TIP 
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FIG  5.  CONTRACTION  OF  ROTOR  WAKE 


FIG.  10.  COMPUTATIONAL  FLOWCHART 
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FIG.  1 1.  AVERAGE  INDUCED  VELOCITY  vs.  INTERATION  NUMBER  - S.58  ROTOR 
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FIG.  14.  VARIATION  OF  BLADE  CIRCULATION  WITH  SPAN  - S.58  HELICOPTER 
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